Enhanced pairing in doped quantum magnets with frustrating hole motion 
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Evidence for strong pairing at arbitrary small J/t is provided in a t-J model on the checkerboard 
lattice for the sign of hopping leading to frustration in hole motion. Destructive quantum interfer- 
ences suppress Nagaoka ferromagnetism when J/t — *■ and reduce drastically coherent hole motion 
in the fluctuating singlet background. It is shown that, by pairing in various orbital symmetry 
channels, holes can benefit from a large gain of kinetic energy. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 



Although the underlying mechanism that leads to 
high critical temperature in the cuprates superconduc- 
tors might be of magnetic origin a detailed scenario is 
still missing. The inter-layer tunnelling mechanism pro- 
posed by Anderson Q is based on the presence of weakly 
coupled CuO layers in this material. The original idea is 
that the pairing mechanism within a given layer (which 
might be rather weak) could be amplified by the Joseph- 
son inter-layer tunnelling of the Cooper pairs. Such a 
scenario is only possible if the single particle motion is 
not allowed along the c-axis perpendicular to the layers, 
i.e. if the coherent single particle tunnelling is blocked by 
some means. This phenomenon could be observed e.g. in 
weakly coupled (along the c-axis) two-dimensional (2D) 
Luttinger liquids 0] (LL) as a result of the orthogonality 
catastrophe 

Despite its beauty and simplicity, this scenario still 
lacks strong experimental support and, on the theoreti- 
cal side, a proof that an analog of the one-dimensional 
LL could be realized in 2D. Nevertheless, I suggest here 
that a closely related mechanism could be realized in 
a doped, strictly 2D, antifcrromagnetic (AF) frustrated 
magnet. In the last years, there has been growing ex- 
perimental and theoretical interests in frustrated mag- 
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FIG. 1: (a) Checkerboard lattice. The plaquette singlets of 
the VBS state are shown schematically in green, (b) Lattice 
isomorphic to the triangular lattice with the same average 
connectivity as (a). 




FIG. 2: (a) Periodic \/32 x \/32 cluster used in this study. 
Note the full Ci v point-group symmetry around the centers 
of the crossed and void plaquettes. (b) Reciprocal space and 
available k-points. Note that the checkerboard pattern of the 
crossed plaquettes leads to a folding w.r.t the dashed line. 

nets which could provide e.g. realizations of new exotic 
gapped spin liquid phases. Transition metal oxides (like 
A2Ti2C-7 titanates [4() with a pyrochlore type structure, 
a three dimensional lattice of corner-sharing tetrahedra, 
exhibit a wide range of interesting physical properties. 
The AF Heisenberg model defined on the checkerboard 
lattice d, @, a 2D analog of the pyrochlore lattice, is 
believed to form a gapped Valence Bond Solid (VBS) Q 
with an ordered arrangement of plaquette singlets on a 
subset of the void plaquettes [f| (see Fig. ^a)) and a 
rather short magnetic correlation length of only a few 
lattice spacings. The GS hence breaks translation sym- 
metry and is two-fold degenerate. 

The recent discovery of superconductivity in 5d tran- 
sition metal pyrochlores Q or in a CoO triangular layer 
based compound |9j suggests that geometric frustration, 
which could be magnetic and/or kinetic, might play a 
key role in the mechanism of unconventional supercon- 
ductivity. In this Letter, I show that pairing appears at 
arbitrary small J/t in the hole-doped checkerboard AF 
magnet for a given sign of the hopping amplitude 
leading to frustration in single hole motion. It is shown 
that hole derealization i.e. gain in kinetic energy can 
play a key role in some mechanisms of unconventional 
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FIG. 3: Single hole dynamics on the checkerboard lattice; 
Spectral functions for t < (a) and t > (b). The total 
weight is a constant (sum-rule) and set to 1; (c) Z-factor at 
the X-point vs J/t; (d) Dispersion width along the k x — k y 
direction in the BZ. For comparaison the case of the square 
lattice is shown. 



pairing, as in the inter-layer mechanism discussed above. 

To describe the low-energy physics of the weakly doped 
magnet the standard t— J model Hamiltonian is used: 
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where the same couplings t and J are used on the ver- 
tical, horizontal and diagonal bonds of the checkerboard 
lattice. Computations are performed by Lanczos exact 
diagonalisation of a 32-site periodic cluster. Although 
tilted by 45° w.r.t. the crystallographic axis, the square 
x V32 cluster (see Fig. exhibits the full C iv 

point group symmetries of the infinite lattice (with 2 sites 
per unit cell) so that all orbital symmetries (s, d x 2_ y 2, 
d xy , etc..) can be considered 0]. In addition, the most 
symmetric k-points are available in the Brillouin zone 
(BZ) as shown in Fig. |2Jb). Note that, since the unit 
cell contains 2 sites, a folding of the Brillouin zone oc- 
curs so that the many-body spectra at k and k + Q (Qo 
is the AF wavevector) are identical (although spectral 
weights are different). An additional folding is expected 
in the thermodynamic limit due to the translation sym- 
metry breaking of the VBS GS. 

Let me first briefly review the properties of a single 
hole doped into the VBS GS of the checkerboard lattice 
(see Fig. |2J) . Numerical computations based on the anal- 
ysis of the hole Green function ^3] (see also Fi gs-Bt a-b)) 
show the existence of quasiparticle (QP) poles [l3j whose 
Z-factor increases with J/\t\ as shown in Fig. |^c). The 
mass of the single hole can be estimated e.g. from the 
dispersion of its QP pole along the k x = k y direction of 
the Brillouin zone (for which the 32-site cluster contains 
5 points). As seen in Ref [13 and shown in Fig. IHId) 
the single hole bandwidth is much smaller than in the 
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FIG. 4: Energies of two holes doped on a square 32-site clus- 
ter (v32x \/32) in various symmetry channels (as indicated). 
The energy reference (Eb ~ 0) is chosen as the bottom of 
the two-hole continuum. Boundstates are characterized by 
Eb < 0. The thick red line corresponds the local (negative) 
slope at J/t — > 0. Thin red arrows show the large (small) 
kinetic (magnetic) energy gain (lost) obtained from the esti- 
mated slope of Eb vs J at J ~ 0.35 (see text). Data obtained 
on a 4 x 4 cluster are also shown using smaller symbols. 



(gapless) square lattice where W ~ 2. 2 J 01 . In the 
latter case, the strong renormalization of the coherent 
motion can be explained e.g. by long-wavelength spin- 
waves scattering |l4|. In contrast, in the checkerboard 
lattice with t > the QP is extremely massive, with 
a very narrow bandwidth, although the VBS host has 
no low energy excitations. This results from a subtle 
interplay between the intrinsic frustrated nature of the 
hopping and the long-ranged plaquette order induced by 
magnetic frustration. Indeed, one expects destructive in- 
terferences 15] between the paths available for the hole 
to hop coherently from one plaquette to the next. How- 
ever, incoherent motion still occurs. In that respect, the 
behavior when J/t — > (for which exact statements can 
be given) is very instructive: while for t < the GS is 
a Nagaoka ferromagnet (of energy — 6|i|), for t > the 
5 = 1/2 GS remains stable down to J = where its en- 
ergy becomes —At (much lower than the fully polarized 
state of energy —2t where the hole is localized on a pla- 
quette). Interestingly enough, the triangular lattice of 
Fig. ^b) exhibits the very same behavior at J/t <C 1. 

As shown below, I argue that holes doped into the 
checkerboard AF magnet pair up very strongly when 
t > 0. The pairing is studied here by considering the 
same 32-site cluster doped with two holes. The two hole 
GS energies are computed for all possible orbital sym- 
metries. The hole-hole binding energies Eb obtained by 
subtracting twice the single hole GS energy are plotted 
in Fig. 01 The results for t > and t < reveal striking 
differences which point toward two completely different 
physical behaviors. The case t < shows a behavior qual- 
itatively (even quantitatively) similar to the case of the 
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well-known square lattice ^(| (for which the sign of t is 
irrelevant): (i) the orbital symmetry of the BS is d x 2_ y 2- 
wave and (ii) the binding energy scales like J (above a 
size dependent critical value) as expected for a purely 
magnetic mechanism. In contrast, the t > case shows a 
completely new behavior: (i) the binding energy appears 
immediately at infinitesimal J; (ii) finite size effects in 
this regime are well controlled: binding is vanishing at 
J = (for all sizes) and the magnitude of the (nega- 
tive) slope characterizing the linear behavior Eb oc J/t 
increases with system size; (iii) the lowest BS has s-wave 
symmetry; (iv) other exotic stable BS appear with non- 
conventional orbital symmetries like d x 2_ y 2-, d xy - or even 
g-wave symmetries. It should be noted that, although 
the s-wave BS wavefunction has the full symmetry of the 
crystal by definition, it could be very anisotropic like the 
superconducting gap in the inter-layer scenario [lj. Im- 
portantly enough, while the BS stability in the thermo- 
dynamic limit at small J/t values (the physical range of 
this parameter) remains controversial in the case t < 
(as for the square lattice 0), the data (and their behav- 
ior with size) for t > provide a very strong evidence for 
the robustness of the BS down to J = 0. Such a remark- 
able phenomenon is due to the nature of the pairing en- 
hancement mechanism which always benefit from a gain 
in kinetic energy. Indeed, using Feynman-Hellmann the- 
orem, the change of magnetic energy in the BS is given 
by Eg ag — J^jf- and hence can be directly estimated 
from the local slope of the function Eb{J) using the sim- 
ple construction shown in Fig. 0] Hence, the difference 
E B ln =E B - E B iag < (oc - J 2 /t when J/t -> 0) should 
be of kinetic origin. It is remarkable that, at small inter- 
mediate J/t couplings (typically ~ 0.2 — 0.6 for 32 sites), 
the stability of the BS benefits only from a large gain of 
kinetic energy while an overall magnetic energy is lost 
in the pairing process. In that regime, the frustrating 
AF magnetic coupling is essential to stabilize the plaque- 
tte VBS structure. Therefore, one expects that such a 
behavior appears for even smaller J/t as singlet plaque- 
tte correlations build up for larger clusters. Then, one 
cannot completely exclude the fact that the J = limit 
becomes singular. 

To investigate further the structure of the two hole BS 
wavefunctions it is of interest to compute the static hole- 
hole correlation Chh( r ), i- e - the probability to find the 
two holes at a given separation r. Results are shown in 
Fig. HJa) as a function of distance. The data for t < 
show a short-distance "depression" in the hole- hole cor- 
relations and a higher probability for the holes to sit at 
distances ~ 3 or 4. As expected, the size of the pair wave- 
function tends to grow as J/\t \ is reduced (similarly to the 
case of the square lattice |l6() as seen also from the be- 
havior of the average hole-hole separation Rhh = \J (r 2 ) 
plotted w.r.t J/\t\ in Fig. 02b). Below a critical J/\t\ 
value, the d-wave BS might then be unstable. In contrast, 
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FIG. 5: Hole-hole density static correlation function Chhijc) 
vs |r| in the 2-hole GS for both signs of t computed on the 32- 
site periodic cluster; (a) Correlations vs distances for J/t = 
0.3 & J/t = 0.6. For t > 0, the data at J = 0.6, close to 
the one shown here for J/t = 0.3, are omitted for clarity; (b) 
Average hole-hole distance Rhh vs J/\t\ for both signs of t. 



for t > 0, the hole-hole correlations in the s-wave BS are 
rather constant with distance showing no tendency for 
the two holes to repel each other (even for small J/t val- 
ues). The fact that a strong signature of an effective 
"attraction" in the quantities plotted in Figs.OJa) & (b) 
is missing is probably directly connected to the origin of 
the pairing which is of kinetic nature and, hence, should 
naturally lead to larger pair sizes. Incidently, one can 
also expect that such pairing states are likely to be less 
sensitive to extended Coulomb repulsion. 

From the numerical data presented above a simple sce- 
nario might be drawn for t > 0. Although its coherent 
motion is very much suppressed, a single hole bears a 
large incoherent motion (see e.g. the large low-energy in- 
coherent weight in Figs.|3fa-b)) and, hence, can melt the 
plaquette VBS in its vicinity. This region, which might 
be fairly extended in space, becomes more favorable for 
a second hole to gain kinetic energy leading to corre- 
lated (or assisted) hopping. It is interesting to point out 
similarities with the case of some frustrated tight binding 
lattices showing both (i) single particle states localised in 
so-called Aharonov-Bohm cages 01 an d (ii) interaction- 
induced delocalized two-particle BS 18]. 

In this Letter, I give evidence for two very different be- 
haviors regarding hole pairing in the t-J model depending 
whether hole motion is frustrated or not, i.e. depending 
on the sign of the hole hopping on a frustrated lattice 
such as the checkerboard lattice. For t < and suffi- 
ciently large J/\t\ values pairing originates from a local 
magnetic effective attraction analogous to the case of the 
non-frustrated square lattice. On the contrary, for t > 0, 
pairing occurs at arbitrary small J/t. The stability of 
the pair BS always benefits from a large gain of kinetic 
energy which might even overcome lost of magnetic en- 
ergy in some regime. In this case, the binding energy 
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FIG. 6: Schematic correlated pair hopping process in the 
checkerboard AF magnet as indicated by red arrows. The 
plaquette units of the VBS are shown in green. Note that the 
pair wavefunction extends on several units (see text). 

reaches a few tenths of t so that high superconducting 
critical temperature would be expected. Interestingly 
enough, pairing occurs in several pairing channels with 
different orbital symmetries. This mechanism bears sim- 
ilarities with the inter-layer mechanism proposed for the 
high temperature cuprate superconductors [lj . Note also 
that, at finite hole density, a paired state with an en- 
hanced kinetic energy gain is very likely to be lower in 
energy than a phase separated state stabilized only by 
a magnetic energy gain. Lastly, I notice that triangular 
and checkerboard lattices behave in a very similar way 
for J ft < 1. 
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